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1. Introduction 

The understanding of the reaUzation of symmetries in QCD from first principles has been 
an important issue for a long time. During the 80's, Vafa and Witten gave arguments against 
spontaneous breaking of parity [jl]] and vector-like global symmetries [Q| in vector-like theories; 
however, these arguments were not as useful as expected. Some years after the publication of 
[|l|], many articles appeared [|3|] calling into question the validity of the paper (see for a recent 
review [Q]). The fact that the issue is still open twenty years after the publication of the first 
paper is an indicative of the complexity of the subject. Regarding the second paper it must 
be remarked that the result is not applicable neither to the Ginsparg-Wilson regularization^ nor to 
one of the most used fermionic regularizations on the lattice for QCD, i.e., Wilson fermions. In 
that case, and as it was shown by Aoki [||, ^, there exists a region of the parameters where parity 
and flavour symmetries are spontaneously broken, for the conditions of the Vafa- Witten theorem 
are not fulfilled in the Wilson regularization. In the end, a theoretical proof of the realization of 
symmetries of QCD is still lacking. 

This is where the p.d.f. formalism can help. The p.d.f. formalism is a powerful tool to analyze 
the symmetries of a theory, widely used in statistical mechanics, and introduced around ten years 
ago in quantum field theories for Grassmann degrees of freedom with success [0]. In this paper, we 
apply the p.d.f. formalism to different regularizations of lattice QCD. The next section is devoted 
to a brief introduction to the p.d.f. formalism. In the second section, we analyze the Aoki phase, to 
find either the existence of a new phase, or a infinite set of sum-rules, the eigenvalues of the Wilson- 
Dirac operator must comply with. The third section applies the same formalism to another system; 
we successfully find, by means of the p.d.f., rigorous proof of parity and vector-like symmetries 
conservation in the Ginsparg-Wilson regularization of lattice QCD at non-zero mass. The last 
section summarizes our conclusions 

2. The p.d.f. formalism 

The usual way to study spontaneous symmetry breaking on the lattice consists in the following 
procedure: An external source, which breaks the analyzed symmetry explicitly, is added. This 
generates a non-zero expectation value of the order parameter for that symmetry. Then we take, 
in this order, first the thermodynamic limit, and finally, the zero external source limit. If the order 
parameter expectation value is non-zero after these two limits, then the symmetry is spontaneously 
broken. Although very popular, the method requires extrapolations to be made. Moreover, in some 
systems, the external source method can not be applied in lattice simulations, for the symmetry 
breaking term may add a potentially problematic sign problem. This is the case of the diquark 
condensate in two colours QCD [||]. It would be desirable to be able to study the fate of the 
symmetries without having to add an external source. 

The p.d.f. formalism enables us to do so. It simply amounts to compute the following quantity 

'The paper ||^ states that vector symmetries are conserved in vector-hke theories, if one is able to find an upper 
bound for the propagator. The paper fails to prove this bound for the Ginsparg-Wilson regularization, as hermiticity 
of the Dirac operator is used during the proof. 
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Pic) = lim ( 5 , 



(2.1) 



with ^{x) the order parameter to be studied. In our case, this is a fermionic bilinear i/a(;c)C?i/a(x), 
where O is a constant matrix. To obtain some computable quantity, we prefer to work with the 



Fourier transform of (2.1) 
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where A is the Dirac operator, and So is the pure gauge action. The expectation values of the 
fermionic bilinear can be computed from (0) easily, taking derivatives of P{q) at ^ = 0, 



d"P{q) 



dq" 



q=0 



dc{icYe'i'P{c] 



■ I [c 



(2.3) 



q=0 



Thus, the moments of the Fourier transform of the distribution function are the expectation values of 
the powers of the observables. For a broken symmetry, the expectation value of the order parameter 
(c) will be zero, for a broken symmetry gives rise to symmetric vacua, and the expectation values 
of the order parameter in those vacua cancel each other. Then, the interesting observables to find 
broken symmetries are (c") with n even. 

3. The Wilson scenario 



Let's apply this machinery to QCD with two degenerated flavours of Wilson fermions [^. 
In this scenario, there exists a phase -the Aoki phase- where parity and flavour are spontaneously 
broken, and this translates into a non-zero value of the fermionic bilinear iyj^t^Y- Surprisingly, 
the expectation value of the bilinear lYTsY equal to zero; this phenomenon was explained in 



10], and a brief hint will be given here: There is a ?7(1) remnant of the original SU (2) flavour 



symmetry, which combines with the original parity operator P to yield zero expectation value of 
iY75¥- other words: We can find a redefinition of parity (a combination of parity and U(l)) 
which remains unbroken. This is the standard picture of the Aoki phase. 

Now we compute these two fermionic bilinears using the p.d.f. formalism: 




(3.1) 



(3.2) 



where Xj the eigenvalues of the 75 A operator at zero external source. The expression (3.1) can be 
easily generalized to any number of flavours. 



We must remark that the p.d.f. can not predict what values will the expressions (3.1) and (3.2) 
take. These depend on the specific properties of the eigenvalues of the chosen discretization of the 
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Dirac operator. In fact, as we will see later, the results change dramatically for the Ginsparg-Wilson 



regularization, even though the expressions (|3JJ) and ( |3.2| ) remain the same 



For ( [3 .11 ) to be zero, as the current picture of the Aoki phase demands, the following cancella- 



tion must happen 




(3.3) 



as we know that the left hand side of the equation must be non-zero, by virtue of (3.2). This 



relationship among the eigenvalues of the 75A operator is non-trivial; in fact, for every even moment 
of the distribution function we obtain a different sum-rule, by enforcing {{Hj/ysY)^") = 
From this point on, we face two different possibilities: 



1. The standard picture of the Aoki phase is right, and these sum-rules must be fulfilled by 
the eigenvalues of 75A. In fact, the p.d.f. can be used to derive easily the sum-rules, which 



remind to those obtained by Leutwyler and Smilga in the continuum in [1 1 1. This point was 



exposed by S. Sharpe in [12] 



2. The current understanding of the Aoki phase is incomplete, for it seems improbable that the 
eigenvalues of 75A comply with {{iij/jsY)^") = for value of n. So there must exist 
a new phase which verifies ((/v/751//)^") 7^ for some value of n. As xPT predicts the 
standard picture for the Aoki phase, the realization of this case would imply that xPT is, in 
some sense, incomplete. An analysis of this point of view was done in [^. 

At this moment, there is no theoretical proof to decide between one of this two realizations. In order 
to distinguish which one occurs, a dynamical fermion simulation in the Aoki phase is mandatory, 
measuring the eigenvalues of the 75A operator, and computing the sum-rules. 

4. The Ginsparg-Wilson scenario 

As we have seen, the original Vafa and Witten theorems fail, due to the existence of exceptional 
configurations, in the later scenario, which, on the other hand, is of paramount importance in lattice 
QCD. So, is there any way we can say something concrete about QCD symmetries?. The answer is 
yes, but we need to choose a 'small eigenvalue free' regularization. It happens that the Ginsparg- 
Wilson^ fermions fulfill this requirement; the p.d.f. will do the rest. 

We denote by D the Ginsparg-Wilson operator; as we know, its eigenvalues Xj are complex, 
and lie in a circumference of radius ^ in the complex plane, whose center is in the real axis, at 
the point ^. Using the standard properties of the Ginsparg-Wilson operator, we can compute the 
eigenvalues of the hermitian operator 75 (D + m), which are 



^^_l±y/{l+am)\lj\^ + m^ Ij ^ ^^^^ 
[ + and/or — [Xj + m) Ay € M 



^Strictly speaking, we will work with a Dirac operator which satisfies {f^'jys} = ays, whereas the Ginsparg- 
Wilson regularization only requires {O^' ,75}= t^Rjs, with R a local operator. Nevertheless, the results can be applied 
to any version of Ginsparg-Wilson fermions, although the calculations may eventually become harder. 
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since when Xj ^ M, the eigenvalues associated to Xj and are paired (±), but if Ay € M, this needn't 
be he case, giving rise to zero modes, and non-vanishing topological charge. As an interesting and 
useful remark, we see that the pL 's are real, and non-zero for a non- vanishing mass m. If we look 
at the expression (pH]), we can see that the modulus of these eigenvalues is bounded from below 
by m (this was remarked in [[l3ll). So we establish that < which in turn implies that, for a 

non-zero mass, the following expectation value^ is zero in the thermodynamic limit 

((.W5T3V^)^) = 2(^^i:-^^<^ v-0. (4.2) 

The summatory is removed, adding a factor equal to the number of eigenvalues 2AV . This result 
states that there is no Aoki phase in lattice QCD with Ginsparg-Wilson fermions. But we do not 
know yet whether the Lagrangian symmetries are spontaneously broken or not. Let's look at the 
next order parameter 



{{iW5V?)=2{-^y-,)-4l{-y-] ). (4.3) 





The first term of the r.h.s. is just equal to ( [4.2[ ), so it must vanish in the thermodynamic limit. The 
second term is minus the square of a real quantity, then it must be negative or zero. The requirement 
(which we will assume) that ixj/ysY be an hermitian operator sets to zero this second term in the 
thermodynamic limit, for the expectation value of the square of an hermitian operator must be 
positive, thence 



: 0. (4.4) 




As both terms in the r.h.s. of (4.3) go to zero as the volume increases, this proves that parity is not 
spontaneously broken in lattice QCD with two flavours of Ginsparg-Wilson fermions^, at least for 
one of the more standard order parameters. In fact, we know that there exists an index theorem for 



Ginsparg-Wilson fermions [|14|], thus we can relate the zero modes of D to the topological charge 
density, 




XT 

Vm^ \ 2 + am 



(4.5) 



Taking into account (p!4]), we deduce that the topological charge density distribution function must 
be a Dirac delta centered in the origin. 



Since we proved that both terms in the r.h.s. of (4.3) must vanish independently in the thermo 



dynamic limit, this result also applies to a single flavoured condensate 



If the expressions ^ and (0), or and Q are compared, one will notice that they are identical, even 
though they correspont to different lattice regularizations. These expressions are regularization independent (see (^). 
But the specific properties of the eigenvalues do depend on the regularization, and make the Wilson and Ginsparg-Wilson 
fermions behave in different ways. 

^The result can be extended to any number of flavours with not much effort. 
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iiiVuYsVu) ) 






(4.6) 



and, by extension, to any linear combination of the single-flavoured condensates xj/jYs^ifj. 

As far as flavour symmetry is concerned, we have proved that ^(/i/Z/sTsV/^)^^ vanish in the 
infinite volume limit, but this is not enough, as this expectation value is forced to be zero because 
of parity conservation. Thus, we would like to investigate the quantity 



2 

y2 




1 



< 



1 



y2 



{Ij + m)' 



[Re{lj)+m]+Im^lj) 



< 



24 
Vm^ 



0. 



(4.7) 



Thus we can affirm that neither parity nor flavour are spontaneously broken in this regularization. 
At zero mass, we cannot establish an upper bound for the observables, hence the argument is not 
valid anymore. The fundamental question is: Why can we bound from above the value of these 
obervables at non-zero fermion mass?. The answer is related to a property of the Ginsparg-Wilson 
operator, that is, {D^^ys} = aRy^, with R a local operator. We can write a similar equation for 
Wilson fermions, where /? is a non-local operator, but in the case of Ginsparg-Wilson fermions, the 
locality of R make the eigenvectors of D look Uke chiral solutions at long distances. So, no quasi- 
chiral, exceptional configurations, at non-zero mass, are allowed, the Aoki phase is completely 
forbidden, and therefore, the symmetries are respected. 

Other interesting results are straightforward from this point on. For instance, we can relate 
the transverse suceptibihty, the topological susceptibility and the chiral condensate. First of all, we 
compute the transverse susceptibility. 



X5 = V{ii\i/Y5Vy 



+ 



4Z7 



^ \ 2 + am 



(4.8) 



Now we write the first summand of the r.h.s. in terms of the chiral condensate, by making use of 
detA = det/sA: 



1 d 9 
— — lndet(Z) + m) 
V dm 




+ 0{a) + 0{ma^ 



(4.9) 



so we amve at 



m 



m 



1-^ 



(4.10) 



where we have dropped the factor (54^) assuming that we are close to the continuum limit. This 
relationship is not new at all, what we are showing here is simply a way to derive it. The interesting 
conclusion comes taking the vanishing mass limit m ^ 0. Then, as the 7] is not a Goldstone boson, 
Xs must remain finite, so 
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lim mxs = ^ lim = y ( W) °^ fWiz- (4.1 1) 

5. Final remarks 

The p.d.f. formalism can be used to cast some light on the old aim of understanding the realiza- 
tion of symmetries of QCD from first principles. Applying the p.d.f. to the Wilson regularization, 
we can explore certain, somewhat overlooked, properties of the Aoki phase. In fact, the p.d.f. 
states that, either the fermionic bilinear iYJsY can take non-zero values in the Aoki phase, ex- 
tending thus the current picture of the phase diagram, or there exists an infinite tower of sum-rules 
the eigenvalues of the Dirac -Wilson operator must comply with. So far, no theoretical argument is 
strong enough to prove one of these scenarios to be right, thus a dynamical fermion simulation is 
mandatory at this point. 

But the most interesting conclusions appear when we apply the p.d.f. formalism to the Ginsparg- 
Wilson regularization. There, we see how parity and vector-like symmetries must be realized for a 
non- vanishing fermion mass. This is a major result that overcomes the difficulties found by 
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